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Abstract. We present X-ray spectral fits to a recently obtained Chandra grating spectrum of 77 Carinae, one 
of the most massive and powerful stars in the Galaxy and which is strongly suspected to be a colliding wind 
binary system. Hydrodynamic models of colliding winds are used to generate synthetic X-ray spectra for a range 
of mass-loss rates and wind velocities. They are then fitted against newly acquired Chandra grating data. We find 
that due to the low velocity of the primary wind (~ 500kms~^ ), most of the observed X-ray emission appears to 
arise from the shocked wind of the companion star. We use the duration of the lightcurve minimum to fix the wind 
momentum ratio at = 0.2. We are then able to obtain a good fit to the data by varying the mass-loss rate of the 
companion and the terminal velocity of its wind. We find that M 2 ~ 10~^ M0yr~^ and Vcx:2 ~ 3000 kms~^ . 
With observationally determined values of ~ 500 — 700 kms~^ for the velocity of the primary wind, our fit implies 
a primary mass-loss rate of A/ 1 ~ 2.5 x 10~* M0 yr~^ . This value is smaller than commonl y infe rred, although 
we note that a lower mass-loss rate can reduce some of the problems noted by Hillier et al. (2001) when a value 



as high as 10"'^ M0 yr~^ is used. The wind parameters of the companion are indicative of a massive star which 
may or may not be evolved. The line strengths appear to show slightly sub-solar abundances, although this needs 
further confirmation. Based on the over-estimation of the X-ray line strengths in our model, and re-interpretation 
of the HST/FOS results, it appears that the homunculus nebula was produced by the primary star. 
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1. Introduction 

The superluminous star tj Carinae (HD 93308, HR 4210) 
continues to be extensively studied over a host of differ- 
ent wavelengths, yet remains intriguingly enigmatic. It is 
amongst the most unstable stars known. In the 1840s, and 
again in the 1890s, it underwent a series of giant outbursts 



5.5 yr period in the variabihty of the He I 10830 A hne. 
Further photometric and radial velocity studies (Damineli 



et al. |l997| , |2000), X- ray observations (Tsuboi et al. |l997 ; 
Corcoran et al. 200C| and references therein), and radio 
data (Duncan et al. 1995, 1999D have supported the 5.5 yr 



(e.g. Viotti 1995) which ejected large masses of material 
into the surrounding medium. HST image s of the result- 
ing bipolar nebula (e.g. Morse et al. 1998 ), known as the 
Homunculus, show it to be amongst the most spectacular 
in our Galaxy. The central object is now largely obscured 
by dust, and the cause of the outbursts and the nature of 
the underlying star (at outburst and today) remain spec- 
ulative. The source continues to show brightness fluctua- 
tions and emission-line variations. Further details of rj Car 
can be found in the review by Davidson & Humphreys 

In recent years evidence for binarity in this system 
has been accumulating. Damineli (1996) first noted a 



period and the binary hypothesis. However, the ground 
based radial velocity curve was not confirmed by higher 
resolution spectra with STIS, indicating that at least the 
time of periastron passage is not well defined by the UV 
and optical spectra (Corcoran et al. 2001a| ). A comparison 
of the abundances from the central object (s) and the com- 
position of the Homunculus nebula has also determined 
th at the re are at least two stars in this system (Lamers et 
al. |l998|) . 
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rj Car is often classified as a luminous blue variable 
(LBV). These are massive stars believed to be in a rapid 
and unstable evolutionary phase in which many solar 
masses of material are ejected into the interstellar medium 
over a relatively short period of time (« 10^ yr). LBVs 
are regarded as a key phase in the evolution of massive 
stars, during which a transition into a Wolf-Rayet star 
occurs (e.g. Langer et al. 1994 ; Maeder & Meynet 2001 ). 
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Due to their rarity and complex nature however, we unfor- 
tunately still have no definitive theory for mass-loss dur- 
ing the LBV stage. The majority of proposed mechanisms 
to drive LBV instabilities, the onset of higher mass-loss 
rates and underlying eruptions, are concerned with the 
importance of radiation pressure within the outer enve- 
lope of the LBV, and for example utilize pulsational in- 
stabilities (e.g. Guzik et al. 199S| ), dynamical instabilities 



1993D , or presuppose Eddington- 



(e.g. Stothers & Chin 

like instabilities. The latter could arise from an enhance- 
ment in opa city a s the star moves to lower temperatures 
(e.g. L amers 1997 ), or fro m the influence of rotation (e.g. 
Langer 1997 ; Zethson et al. 1999). Alternatively, the possi- 
bility that binarity plays a fundamental role in explaining 
observed LBV outburst properties has also been consid- 
ered (Gallagher 1989 ), though most LBVs are not known 
binaries. Clearly, determining the wind and stellar prop- 
erties of LBV stars is paramount (see, for example, the 
discussions in Leitherer et al. 1994 and Nota et al. 1996 ). 

An important question is the degree to which binarity 
influences the properties of LBVs (i.e. do LBVs in bina- 
ries evolve differently than single LBVs?). So while the 
presence of a companion can be exploited to help mea- 
sure the mass of such stars, we must bear in mind that 
binary LBVs and single LBVs may be quite distinct ob- 
jects. Therefore, in order to use rj Car to understand some 
of the defining LBV characteristics such as their extremely 
high mass-loss rates, we first need to determine beyond all 
doubt that rj Car is in fact a binary, and then to determine 
the influence of the companion on the system. 

Investigations over the last few years have already 
helped to form a basic picture of ij Car. The orbital 
parameters, although uncertain, indicate the presence of 
an early-type companion star, which will also have a pow- 
erful stellar wind. In such binaries, a region of hot shocked 
gas with temperatures in excess of 10 million K is cr eated 
where the stell ar win ds collide (Prilutskii & Usov 1976| ; 
Cherepashchuk 1976 ). The wind- wind coUision (WWC) 
region is expected to contribute to the observed emission 
from this system, particularly at X-ray and radio wave- 
lengths. Previous X-ray observations revealed extended 
soft emission from the nebula and strong, hard, highly ab- 
sorbed, and variable emission closest to the star (Corcoran 
et al. 1995; Weis et al. 2001), in contrast to the emis- 
sion characteristics from single stars, which are typically 
softer, much less absorbed, substantially weaker and rel- 
atively constant in intensity. Since 1996 Feb t] Car has 
been continuously monitored by RXTE in the 2-10 keV 
band (e.g. Corcoran et al. 2001a). The lightcurve (Fig. |l]) 
contains remarkable detail showing a slow, almost linear, 
rise to maximum over a period of « 1 yr, followed by a 
rapid drop to approximately 1 /6 of the peak intensity for 
« 3 months, an almost as sharp rise to approximately 1/2 
of the peak intensity level, and then almost constant in- 
tensity for « 3/4 of the proposed 5.5 yr orbital period. 
The drop to minimum was successfully predi cted from 
numerical models of the WWC (Pittard et al. |l998|) be- 
fore being actually observed. 



Small scale quasi-periodic outbursts in the X-ray 
lightcurve have also been detected (Corcoran et al. 1997). 
Estimates of changes in the timescale between succes- 
sive flares as a function of phase were made by Davidson 
et al. (1998) for a variety of assumed orbital elements. 
RXTE X-ray observations obtained after the X-ray min- 
imum seem to show a lengthening of the flare timescale 



(Ishibashi et al. 1999), which indirectly support the binary 
model. The latest published X-ray observation of rj Car 
is of a high resolution grating spectrum taken with 
the Chandra X-ray observatory (Corcoran et al. 2001b). 



Preliminary analysis has revealed the presence of strong 
forbidden line emission, which suggests that the density 
of the hot gas, < lO^^'cm"'^. This can be contrasted 
with the newly published X-ray gratin g spectra of the sin- 
gle stars 9^ Or i C (Schulz et al. |2000| ), ( On (Waldron & 
Cassineni |200ll ), and C Pup (Kahn et al. ^00l| ), aU of which 
have weak forbidden lines (indicative of either high den- 
sities or high UV flux near the line forming region). This 
lends further support to a wind-wind collision model, al- 
though it is possible that some of the forbidden emission 
may be related to the surrounding nebula. 

If 77 Car is in fact a binary system, the orbital elements 
and stellar parameters are not yet tightly constrained 
by current observations. Ground-based observations (for 
which good phase coverage exist) are hampered by poor 
spatial resolution and thus suffer contamination from 
strong nebular emission. High-spatial resolution spectra 
have been obtained by STIS but phase coverage is cur- 
rently very limited. HST STIS observations at two differ- 



ent phases of the 5.52 year cycle (Davidson et al. 200C) did 
not confirm the predicted variations in the radial velocity 
of the emission lines based o n the ground-based radial 
velocity curve (Damineli et al. 2000 ). If 77 Car is a binary, 
it is vitally important to determine the stellar parameters 
of the companion so that the effect of the companion on 
observations can be understood and the correct stellar pa- 
rameters of the primary can be derived. 



1.1. X-ray Emission from Colliding Winds 

The wealth of information contained in X-ray spectra 
of colliding wind binaries (e.g. the density, temperature, 
velocity, abundance, and distribution of the shocked gas 
in the wind collision region) has been a strong motivat- 
ing force for observers and theorists alike in this field. 
Since the hot plasma in most colliding wind shocks is op- 
tically thin and coUisionally ionized, and is generally as- 
sumed to be in collisional equilibrium, Raymond-Smith 
(Raymond & Smith |1977| ) or MEKAL (Mewe et al. [1995D 
spectral models are normally fitted to such data (e.g. 
Zh ekov fc Skinner 2000| ; Rauw et al. ^000 ; Corcoran et 
al. 2001b). However, the multi-temperature, multi-density 
nature of the WWC region means that at best simple fits 
with one- or two-temperature Raymond-Smith models can 
only characterize the broad properties of the emission. In 
this way one can estimate an 'average' temperature of the 
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Fig. 1. Lightcurve of 77 Car observed with the RXTE satellite and phased to the 5.5 yr orbital period. Plotted are counts 
detected in layer 1 of the second proportional counter unit (PCU2) and a predicted lightcurve (Pittard et al. 1998) 
from a numerical model of the wind-wind collision. The two agree well, particularly the duration of the minimum. The 
rise from minimum is not in good agreement, but this is thought to be due to the limitations of modelling the wind 
collision in 2D. The rapid change in position angle of the stars through periastron passage skews the shock cone which 
causes the line of sight in fully 3D mo dels t o remain in the denser wind of the primary until later phases, increasing 
the absorption at these times (Pittard 2000) 



shocked gas, and an 'average' absorbing column, but lit- 
tle is learned of the underlying stellar wind parameters. 
At worst the application of one- or two-temperature mod- 
els to what is inherently multi-temperature emission can 
lead to spurious values of some of the fit parameters, e.g. 
abundances (cf. Strickland & Stevens 1998). 



Complex numerical hydrodynamical models have often 
been applied to gain insight into colliding wind systems 
(e.g. Stevens et al. [l992|, Pittard et al.|l998|). However, 



while undoubtedly useful, their interpretation can be 
difficult, and to date there have been only two published 
papers where observed spectra are directly fitted with syn- 
thetic spectra from such models. In the pioneering work 
of Stevens et al. ( 1996| ), medium resolution ASCA spectra 
of the Wolf-Rayet binary 7^ Velorum were fitted against a 
grid of synthetic spectra. In this fashion they were able to 
obtain direct estimates of the mass-loss rates and termi- 
nal velocities of the individual stellar winds. As mass-loss 
rates obtained from measures of radio flux or spectral line 
fits depend on a variety of untested assumptions, the im- 
portance of a new independent method to complement 
estimates from free-free radio or sub-mm observations, or 
from Ha or UV spectral line fitting, cannot be stressed 
enough. Rates derived by Stevens et al. (1996) with this 
new method were significantly lower than the commonly 
accepted estimates for 7^ Velorum based on radio obser- 
vations, but an indication of the future benefits of this 
method was realized when both sets of estimates were later 
brought into agreement following a surprisingly large re- 



duction in the distance to this star from Hipparcos dataQ 
We note that this method can also provide insights into 
the values of parameters which are otherwise virtually im- 
possible to estimate, such as the mass-loss rate of the com- 
panion, M2, or the characteristic ratio of the prc-shock 
electron and ion temperature (Zhekov & Skinner 2000). 
The quality of recently available X-ray grating spectra 
now gives us access to important X-ray emission line di- 
agnostics which should severely constrain models of the 
X-ray emission distribution. This means that stellar wind 
parameters can in principle be reliably estimated from 
analysis of X-ray grating spectra of colliding wind bina- 
ries. 

In addition to testing the binary hypothesis, the 



Chandra grating spectrum of Car (Corcoran et al. 2001b) 
provides the ideal opportunity to test the method de- 
veloped by Stevens et al. (1996) against a spectrum of 



much higher spectral resolution, and to pin down impor- 
tant physical parameters of the system. In this paper we 
therefore fit the X-ray grating spectrum using a grid of 
colliding wind emission models to i) test the binary hy- 
pothesis, and ii) to attempt to obtain accurate estimates 
of the wind parameters of each star. The fact that we 
are able to obtain good fits, with sensible model param- 
eters, gives us further confidence in the binary hypoth- 
esis. We also find that unlike the UV and optical work 



^ The thermal radio flux, oc A/"/^!)"^ (where D is the 
distance to the source) whereas the X-ray flux from an adia- 
batic wind collision is oc M ^D~^. Therefore M radio oc D^^^ 
whereas Mxray oc D. If 
creases faster than Mxr, 
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where fits to the primary are made difficult by signifi- 
cant contamination from the companion star, instead the 
X-ray emission arises from the shocked wind of the com- 
panion and suffers essentiaUy zero contamination from the 
wind of the primary. Therefore the X-ray data uniquely 
samples parameters of the companion, in contrast to the 
optical analysis which probes the nature of the primary. 
In this sense our analysis is entirely complementary to 
the complex fits to the UV and optical HST spectrum 



of T] Car by Hillier et al. (2001). Our analysis also pro- 
vides us with a new estimate of the mass-loss rate of the 
primary star. For details of the Chandra observation and 
an initial analy sis of the data the reader is referred to 
Corcoran et al. ( ^001b| ). Here we only note that there is no 
significant contamination of the dispersed spectrum from 
any spatially resolved emission (i.e. in the Homunculus). 
In Section g we discuss the creation and variation of the 
model spectral grid; in Section^ we describe the fit results; 
and in Section ^ we summarize and conclude. 

2. The Synthetic Colliding Wind Spectra 

The method applied by Stevens et al. ( 1996|) as discussed 
in Section [| is as follows. We first calculated a whole 
host of hydrodynamical models with different stellar wind 
parameters. We then generated a synthetic X-ray spec- 
trum from each model. Finally, we fit the grid of synthetic 
spectra to the actual data. Best-fit values are then ob- 
tained for each parameter of interest (e.g. Mi, M 2, Wooi , 
V0C2 etc. ). Specific points in relation to this method are 
highlighted below. 

2.1. Hydrodynamical Models off] Carinae 

The colliding wind models were calculated using VH-1 
(Blondin et al. 1990 ), a Lagrangian- remap version of the 
third-order accurate Piecewise Parabolic Method (PPM; 
Colella & Woodward 1984 ). The stellar winds are mod- 
elled as ideal gases with adiabatic index 7 = 5/3. 
Radiative cooling is included via the method of opera- 
tor splitting, and is calculated as an optically thin plasma 
in ionization equilibrium. The cooling curve for the tem- 
perature range 4.0 < log T < 9.0 was generated using 
the Raymond-Smith plasma code. Despite some evidence 
to the contrary (see Corcoran et al. 2001b), we have used 
solar abundances throughout this work, since it keeps the 
problem as simple as possible in this first investigation. 
We aim to model non-solar abundances in a later paper. 

The simulations were calculated assuming cylindri- 
cal symmetry of the wind collision zone - orbital effects 
are negligible at the phase of the Chandra observation 
{(j) = 0.60). Grid sizes spanned a range from 300 x 300 to 
550 X 450 cells. Each grid cell was square and of constant 
size on an individual grid. The linear dimension of each 
cell was either 2.67 x 10^^ or 4.0 x 10^^ cm. The maximum 
distance from the axis of symmetry was 1.2 x 10^^ cm 
in all cases. At the phase of the observation, the orbital 
separation of the stars using the ephemeris of Corcoran et 



al. ( |2001aD is « 4x lO" cm. Thus there is either 100 or 150 
grid cells between the stars. The large separation means 
that the winds are likely to collide at very near their termi- 
nal velocities. We therefore assume that we can treat the 
winds as being instantaneously accelerated to their termi- 
nal velocities, and do not consider any radiative driving ef- 
fects. In comparison with the work on 7^ Velorum (Stevens 
et al. 1996), assuming terminal velocity winds and negli- 



gible binary rotation is more valid for 77 Car. We further 
assume that the winds are spherically symmetric. 

To compute a grid of synthetic spectra we initially 



varied four parameters (Mi, M 2, 



the 



hydrodynamic models. During our investigation we also 
included the separation of the stars, D, as an additional 
free parameter. From fits at this point it was found that 
_D w 4 x 10^^ cm with a small uncertainty, in good agree- 
ment with the expected value from our current under- 
standing of the orbit. We therefore fixed it at this value 
for the rest of our analysis. However, we found that there 
were large uncertainties on the values of Mi and Wooi- 
With the benefit of hind-sight this is not too surpris- 
ing given the known slow speed of the primary wind 



(Wo 



500 — 700 km s ^ ) , which means that the shocked 



primary wind is not a strong source of X-ray emission at 
hard energies. Therefore for our final grid we fixed the 
terminal velocity of the primary star at 500 kms~^ and 
adjusted Al 1 to obtain a desired value for the wind mo- 
mentum ratio, T] = M 2 i'oo2/(-^i ''^ooi), with M2 and 
V002 as free parameters. 

The range of the free parameters used in our final grid 
is given in Table ^ To restrict the number of models to a 
manageable number the parameter steps are fairly coarse. 
In future papers we will use a much finer grid. The range 
in the value of 77 corresponds to either the wind of the 
primary (77 — 0.1) or the secondary (77 — 5) dominating, 
or the winds having equal momentum fluxes (jj = 1). The 
distance of the stagnation point from the centre of the 
primary star is given by 



1 + 



(1) 



and ranges from 1.2 x 10^^ cm to 3.0 x lO^'' cm. In the 
most extreme cases, this is still of order ten to a hundred 
times the radius of the stars, and justifies our use of termi- 
nal velocity winds. The kinetic power of the primary wind 
ranged from 4.7 x 10'^'* ergs~^ to 7.9 x lO'^^ ergs""'^ , and 
from 7.1 X 10^^ ergs^^ to 7.9 x 10^* ergs~^ for the sec- 
ondary wind. The combined kinetic power of the winds 
ranges from 7.6 x 10'^'^ ergs^^ to 1.6 x 10"^^ ergs"^ . The 
half-opening angle of the contact discontinuity measured 
from the line between the secondary star and the shock 
apex ranges from « 50 - 120° (cf. Eichler & Usov |1993D . 

The effect of radiative cooling can be quantified by 
the parameter x, the ratio of the cooling timescale to the 
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Table 1. Our final grid of stellar wind parameters for 
the hydrodynamical models used to generate synthetic 
spectra. The models fix the terminal velocity of the pri- 
mary wind at 500 kms^^ and the stellar separation at 
D = 4.0 X 10^'' cm. The mass-loss rate of the primary star 
is given by Af i = M 2 Wooa / (500 rj) . 



Parameter 


1 


2 


3 


V 


0.1 


1.0 


5.0 


M2 (Moyr-i) 


10"'' 


10"^ 


10-^ 


Uooa ( km s ^ ) 


1500 


3000 


5000 



dynamical timescale of the system. For shocked gas near 
the local minimum in the cooling curve at T « 2 x 10^ K, 



X 



M ^7' 



(2) 



where Wg is the wind velocity in units of 1000 kms"^ , 
di2 is the distance to the contact discontinuity in units 
of 10^^ cm, and M _7 is the m ass-los s rate in units of 
10"'^ Mgyr^i (cf. Stevens et al. |l992|) . This equation is 
valid for post-shock temperatures in the range 10^ ^ T < 
lOS-2 K (680/^ < < 2700/Vm kms"! where ^mn is 
the average particle mass in grammes) . For material with 
solar abundance (/i = 0.6), this corresponds to 0.9 < ws ^ 
3.5. For winds with slower velocities, the shocked gas lies 
on the negative slope of the cooling curve. In this case the 
dependence of x on the velocity of the wind is steeper: 



X ■ 



M_7 ■ 



(3) 



In our grid of simulations, the shocked wind of the primary 
star is almost always strongly radiative (xi,min << 1-0) 
while the secondary's is never strongly radiative (x2,min — 
2.0). However, for some choices of the 3 parameters in 
Table |l|, the primary's wind can approach the point where 
it starts to become adiabatic (xi,max ~ 1-0). In Fig. ^ we 
show two examples of the hydrodynamic calculations used 
to calculate synthetic spectra. Due to the difference in ter- 
minal velocity of the two winds, the velocity shear at the 
wind interface always generates Kelvin-Helmholtz insta- 
bilities in our models. The rapid cooling of at least one 
of the two winds means that thin-shell instabilities also 
always occur in our models. The colliding winds region is 
most unstable when both shocked winds rapidly cool (cf. 
Stevens et al. [19921 ). 

The spectrum from each hydrodynamic model was 
averaged over 3 'snapshots' each spaced by about 70 d. 
This is of order the wind dynamical timescale and serves 
to approximate a time-averaged spectrum. It is particu- 
larly important to adopt this approach for models where 
the wind collision region is very unstable since in these 
cases the flux can vary by greater than ±50% from one 
sna pshot to the next. As already noted by Corcoran et 
al. (2001b) no variability is seen during the « 90 ksec ex- 
posure of the Chandra grating spectrum. This is to be 



expected, however, since the dynamical timescale is much 
longer. Two examples of 'time-averaged' synthetic spectra 
are shown in Fig. ||, next to density plots of the corre- 
sponding hydrodynamic calculation. 

The synthetic spectra were calculated over the energy 
range 1.26-10 keV (10-1.26 A respectively), and have a 
resolution of 0.005 A. The actual grating spectrum shows 
significant absorption below 1.5 keV and contains es- 
sentially no useful information below our lower limit of 
1.26 keV. The spectral resolution of our synthetic spectra 
is approximately twice as high as the grating data that we 
model. Use of the Raymond-Smith plasma code implicitly 
assumes that the plasma is in thermal equilibrium. This is 
generally true for colliding wind binary systems (see Luo 
et al. 1990) and is a good assumption for 77 Car. 



2.2. Variations with i], M 2, snd v, 



002 



In Fig. H and the following subsections we reveal how the 
shape and normalization of the calculated spectra depend 
on the various free parameters of the grid. 



2.2.1. Variation with r] 

Due to the low preshock velocity of the primary wind, the 
X-ray emission from the WWC is almost entirely from 
the shocked secondary wind. A greater fraction of the sec- 
ondary wind passes into the wind collision zone as the 
relative momentum of the primary wind is increased, so 
we expect the luminosity to increase as 77 decreases, as in- 
deed the top panel of Fig. ^ shows. An important point, 
however, is that the shape of the spectra do not seem de- 
pendent on the value of ry, contrary to our expectations 
before starting this study. 



2.2.2. Variation with M2 

As we increase the mass-loss rate of the secondary star in 
our models, the kinetic power and the density of its wind 
also increase. Therefore, for a given wind momentum ratio, 
the rate of conversion of kinetic energy to thermal energy 
increases, as does the density of the postshock gas. Thus 
there is more available energy to radiate and a greater 
efficiency of radiation. Hence we should expect the X-ray 
luminosity to scale strongly with increasing M 2- This is 
indeed what we see in the middle panel of Fig. where an 
order of magnitude increase in M 2 leads to an almost two 
order of magnitude increase in X-ray luminosity, indicat- 
ing that the calculations in this panel are in the adiabatic 
regime {Lx oc M^; see Stevens et al. 1992). However, there 
is no evidence for a softening of the spectrum with increas- 
ing M 2- This is again contrary to our initial expectations 
(note, however, the discussion in the next section). 
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Fig. 2. Hydrodynamic simulations and theoretical spectra of the colliding winds in 77 Car. The top panels show 
the results from a model where the momentum of each wind is equal (?/ = 1), whereas in the bottom panels the 
primary wind has a momentum 10 times greater than the secondary wind (77 = 0.1). In both cases a density plot 
(g cm~^) from the hydrodynamical grid is shown on the left and the resulting intrinsic X-ray spectrum shown on 
the right (no absorption). On each hydrodynamical grid the primary star is located at (0,0) and the secondary star 
is at (4.0 X lO^^'jO). The shape and structure of the shocked region varies with the momentum balance between the 
winds and the respective value of x for each shocked wind. The full parameters used in these models were: upper 
panels - Mi = 10"^ Moyr'^ M 2 = 10^^ Mgyr'^, Vaoi = 500 kms"S ^002 = 5000 kms~^ ; lower panels - 



M 1 = 3 X 10"'' Mq yr"^ , M 2 = 10"^ M© yr^ 



500 kms 



1500 kms . The luminosity and shape 



of the X-ray spectrum are a direct consequence of these parameters. 



2.2.3. Variation with V002 

If the preshock velocity of a wind is increased, the post- 
shock temperature increases (T cx w^). Therefore we ex- 
pect the model spectra to harden as the preshock velocity 
of the secondary wind is increased. Because the kinetic 
power of the wind increases we also expect the luminosity 
to increase. The variation of the spectrum with V002 is 
shown in the bottom panel of Fig. ^ Here we find that our 
expectations are only partially met. The overall luminosity 
does increase with V002 > and the spectrum does harden be- 
tween V002 = 1500 — 3000 kms~^ , but it softens between 
1^002 — 3000 — 5000 kms^^ . An investigation into this 
surprising finding offers the following explanation. 

For fixed 77 and M 2 , and variable V002 j the density of 



the wind varies as p oc M / V002 ■ If ^'002 is increased p de- 
creases, which reduces the radiative efficiency of the hot 



shocked gas near the stagnation point (E cx n^), which 
acts to slightly suppress the hard flux. A plot of the aver- 
age density against T for tioo2 = 1500, 3000, 5000 kms~^ 
also reveals that for 10^'^ < T < 10'''' the average 
density is higher when t;oo2 = 5000 kms""' than when 
V0D2 = 3000 kms~^ , which tends to keep up the soft 
emission relative to the hard emission. This appears to 
be because there is some mixing at the interface be- 
tween the primary and secondary wind, and the shocked 
primary wind is denser when V002 — 5000 kms^^ 
(Ml = 3 X 10-^ 6 X 10-^ 10-3 M0yr-i for V002 = 
1500, 3000, 5000 kms"! ). The mixed gas tends to be at 
these intermediate temperatures. 

The overall net effect is that the emission from the hot 
gas is suppressed relative to the cooler gas, which softens 
the spectrum from the model with Vao2 — 5000 kms~^ 




Fig. 3. Variation in the spectral shape and normahzation of the theoretical colliding winds spectra in rj Carinae. 
The top panel shows the variation with 77 {M 2 = 10^^ Moyr"^ , V002 = 3000 kms^^ ). The middle panel shows 
the variation with M2 {1] = 5.0, V0C2 = 1500 kms^^). The bottom panel shows the variation with Vco2 iv = 0.1, 
M2 = 10"^ Meyr-i). 



relative to the model with V002 = 3000 km s ^ . These in- 
ferences are supported by spectra calculated with U002 — 
5000 kms-i and r] = 0.1, and variable M2 (Fig. |). This 
shows that the spectrum hardens as M 2 increases from 
10"^ Mgyr-i to 10"^ Mgyr"! (the spectral index a, 
where oc v", increases by about 0.1). This is what 
we expect from our prior inferences since a higher M 2 
helps to make postshock gas with a preshock velocity of 
V002 = 5000 kms~^ somewhat more efficient at radiating. 



In summary, we find that varying M 2 can actually have 
a small effect on the shape of the spectrum in some parts 
of parameter space, as well as the more obvious and much 
greater effect on the luminosity. The power of M 2 on af- 
fecting the spectral shape is, however, much less than that 
of V002 1 which is the primary influence. 
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Fig. 4. Variation in the spectral shape and normahzation 
for a very fast secondary wind (woo2 = 5000 kms~^ ) and 
a fixed wind momentum ratio [rj — 0.1). As M 2 increases 
the spectrum becomes shghtly harder. 



2.2.4. Summary of Spectral Dependence on rj, M 2, 

and V002 

For the models investigated in this paper, we find that 
the spectral shape is i) insensitive to the value of rj, ii) 
can harden slightly as M 2 increases (in other parts of 
parameter space, e.g. strongly radiative shocks, it may 
soften with increasing M2), and iii) generally hardens 
with increasing Vo^^ i but can sometimes soften. The 
luminosity increases with i) decreasing rj (since 77 Car is 
unusual in the sense that there is no discernible contri- 
bution to the X-rays from the shocked primary wind, bar 
mixing), ii) increasing AI 2, and iii) increasing Woo2- 

2.3. The Spectral Grid 

The synthetic spectra were input into a FITS-style 'Table' 
model for use with t he XS PEC data-analysis package, as 
specified by Arnaud (1999). The fitting procedure within 
XSPEC is able to interpolate between the parameter val- 
ues on the spectral grid so that best-fit values of M2 
etc. are obtained from fitting the Chandra data. In the 
2 — 10 keV band, the intrinsic luminosity on our grid 
ranges from 1.8 x 10^^ ergs~^ to 1.6 x 10^^ ergs^^. 



which brackets the observationally determined value of 
4 X lO^-* ergs-i (Corcoran et aL ^OOlhl ). 



3. Fits to the Chandra grating spectrum 

We compared the model intrinsic colliding wind spectra 
with the observed grating spectrum allowing the mass- 
loss rate and wind terminal velocity of the companion to 
vary. We in addition included absorption from a cool over- 



lying medium to simulate the combined absorption of the 
circumstellar (nebular) material in the Homunculus, the 
ISM, and the wind from the companion (all of which may 
contribute to the X-ray absorption at some level). Since 
the stellar wind absorption depends on the companion's 
mass-loss rate (at the phase of the Chandra observation) 
for consistency we should model this component of the 
absorption directly, but this is beyond the scope of the 
present paper. Because of the low number of counts per 
spectral bin in the high-resolution grating spectrum, and 
because background does not contribute significantly to 
the observed spectrum at energies above 2 k eV, w e fit the 
gross spectrum using the C-statistic (Cash 1979 ), which 
is appropriate for Poisson-distributed data. 

The fact that the shapes of our model spectra are not 
dependent on the value of 77 introduces a degeneracy into 
our grid as both 77 and AI 2 primarily influence the nor- 
malization of the spectra. This means that various com- 
binations of 77 and AI 2 can provide similar quality fits to 
the data. However, it is possible to obtain an estimate for 
rj from the duration of the lightcurve minimum (sa 100 d). 
To do this we have constructed a simple model of the ob- 
served X-ray emission. We assume that the intrinsic flux 
varies as l/D, and that the absorption from the shock 
apex along the line of sight is negligible when viewing 
through the less dense secondary wind, but total when 
viewing through the very dense primary wi nd. We use 
the latest orbital parameters (Corcoran et al. 2001a[) and 
assume that the skew angle of the shock cone from the or- 
bital velocity, t^orb, is (5 ~ arctan(7;oib/i^ooi)- Fig. ^ shows 
the results, scaled so that L^, = 1 at periastron. The skew 
of the shock breaks the symmetry of the observed emis- 
sion so that the post-minimum flux is lower than the pre- 
minimum flux, as observed. The duration of the minimum 
decreases with increasing rj (for 77 = 0.1,0.2,0.3 the du- 
ration is 133,92,71 d), and is best matched hy rj O.2Q. 
Therefore we adopt this value for the rest of our analysis. 

We attempted to fit the grid of synthetic colliding 
wind spectra to the observed MEG -f 1 order spectrum. 
We also constrained the spectrum normalization to lie be- 
tween 0.8 and 1.2, which fixes the distance to the star 
between 2300 and 1900 pc, close to the canonical distance 
of 2100 pc. In Fig. ^ we show the best fit interpolated 
spectrum from our grid in the range 1.5 — 7.5 keV. Fig. |^ 
shows a close-up of the fit to the Si and S emission lines 
between 1.8 — 2.8 keV. Both figures demonstrate the ex- 
cellent agreement found between the data and the mod- 
els and highlight the progress made in both the observa- 
tional and theoretical study of colliding winds over the 



^ Incidently, the lightcurve minimum computed with the 2D 
hydrodynamical model in Pittard et al. 199J (redisplayed in 



Fig. |l| at the beginning of this paper) also matches the dura- 
tion of minimum very well (although not the post-minimum 
flux because the skew of the shock could not be computed 
with a 2D calculation) . Although different mass-loss rates and 
a slightly less eccentric orbit were used, the wind momen- 
tum ratio adopted was also 77 — 0.2, which provides a certain 
robustness for this value. 
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Fig. 5. The duration of the X-ray minimum predicted 
from simple models. The wind momentum ratio in the 
three models shown is -q — 0.2,1.0,5.0 (solid, dots, 
dashes), which gives a duration of 92, 33, 16 d respectively. 
The observed duration is « 100 d, so the data is best ap- 
proximated by rj = 0.2. 

Table 2. The best fit parameters from the fit of the spec- 
tral grid to the grating data. We fix the terminal velocity 
of the primary wind at 500 kms~^ , the stellar separation 
at D — 4.0 X 10^"* cm, and the wind momentum ratio (77) 
at 0.2, and allow M 2, ^0021 -^h and the normalization to 
vary. Solar abundances were assumed, and a cold external 
absorbing column was applied. The resultant normaliza- 
tion is very close to 1.0, which indicates that the model 
and the values of D and rj are sensible. Alternatively, it 
implies that our assumed distance to rj Car is too far by a 
factor of 1.077, if we imagine this as the sole source of dis- 
crepancy. This would revise the source distance down to 
2.1/1.077 = 1.95 kpc, which is well inside of observational 
uncertainties. The implied mass-loss rate of the primary 
star from our fit is M 1 = 2.5 x 10"** Mq yr~^ . 



Parameter 



Value 



V 

M2 (Meyr-i) 
V002 (kms"-^ ) 
A^H (cm-2) 
Normalization 



0.2 (fixed) 



3+0.09 
^-0.08 

7.68«:- 
1.16 



X 10^ 



last 10 years. The continuum shape and level from the 
models matches closely that of the data. All strong lines 
(and most weak ones) which appear in the observed spec- 
trum are matched also in the best fit. This means that the 
temperature distribution seen in the grating spectrum is 
matched by the model. However most of the strong lines 
are significantly overpredicted by the model fit. This dis- 
crepancy is in the opposite sense to the fit results from 
simple one- and two-temperature Raymond-Smith mod- 
els, and is perhaps revealing that our assumption of solar 
abundances needs to be modified. The best fit parameters 
are summarized in Table 3. 



4. Conclusions 

In this paper our aim has been to test the binary hypoth- 
esis of ?7 Car by directly fitting its X-ray spectrum using a 
grid of spectra calculated from hydrodynamical models of 
the wind-wind collision. While our analysis does not prove 
that it is a binary, we find that the colliding wind emission 
model naturally provides for the range of ionization seen 
in the emission line grating spectrum for reasonable values 
of the wind parameters. We have not shown that it is in- 
consistent with emission from a single star, but have noted 
that it is unlike any of the other single stars observed so 
far at high energies and dispersion. 

The technique applied in this paper has only been 
demonstrated once before and is the first time that it has 
been used with a high quality grating spectrum. Due to 
the low velocity of the primary wind (« 500 km s^ ^ ) , most 
of the observed X-ray emission arises from the shocked 
wind of the companion star. We find it difficult therefore 
to fit both M 1 and Vooi as free parameters. However, the 
duration of the observed X-ray minimum can be used to 
estimate the wind momentum ratio of the stars, rj. With 
r] fixed at 0.2, and M 2, and Wooa as free parameters, we 
are able to obtain a good fit to the data. 

We find that the mass-loss rate of the companion is 



M 



and the terminal velocity of its 

1 



M2 ~ 10 

wind is W002 ~ 3000 kms~^ . These values suggest that 
the companion is probably an Of supergiant (0-stars 
with similar wind parameters - e.g. HD 15570 (04If), 
HD 93129A (03If), HD 93250 (03Vf), HD 151804 (08If), 
and C yg 0B2 #9 (05If) - are fisted in Howarth & Prinja 
1989 ), or is possibly a WR star. The velocity of the pri- 
mary wind has been determined to lie in the range 500 — 
700 kms~^ (e.g. Hillier et al. 2001). Hence our fit implies 
a primary mass-loss rate of M 1 « 2.5 x 10"'' Mq yr~^ . 
From the uncertainty in the value of 77, and the interpo- 
lation on our grid, we estimate the uncertainty on our 
derived value for M 1 as approximately a factor of 2. 

Our best-fit estimate of M 1 is smalle r tha n typi- 
cally infer red (cf. Davidson & Humphreys 1997; Hillier 
et al. 2001) However, we note that a lower mass-loss rate 
for the primary star can reduce some of the problems 
noted by Hillier et al. (2001) who fitted a value as high 
as 10^"^ M© yr~^ . In particular, the models of Hillier et 
al. (2001) suffered from absorption components that were 
too strong and electron-scattering wings which were over- 
estimated. Both indicate that their chosen mass-loss rate 
is too high. Paradoxically both the Ha and H/3 emission 
lines were weaker than observed, indicating that their 
mass-loss rate is too low. However, it is well known that 
the wind collision zone can be a strong source of emission 
lines (e.g. HD 5980, Moffat et al. [l998D , which would re- 
solve this problem. Their inferred minimum column den- 
sity is also larger than the observed X-ray value, again 
implying an overestimate of M 1 . To increase the mass- 
loss rate of t he pr imary towards the value estimated by 
Hillier et al. (2001), we require either a reduction in the 
wind velocity to w 100 kms~^ , which in turn is in con- 




^1 2 3 4 5 6 7 8 

channel energy (keV) 

Fig. 6. Best fit to the Chandra grating spectrum. The fit parameters were M2 — 0.98 x fO^^ Moyr^^, U002 = 
3000 kms^^ , T] ~ 0.2, A'^h = 7.7 x 10^^ cm~^, and the overall normalization was f .16. If we take Wooi = 600 kms~^ , 
this implies that Mi — 2.5 x 10^^ M0 yr~^ . As seen from the figure, the continmmi level is very well fitted by the 
interpolated model. 



flict with previous estimates, or a reduction in the wind 
momentum ratio to 77 ^ 0.1, which is in conflict with the 
observed duration of the X-ray minimum. Finally, it is 
worth noting that our results indicate a value for M 1 
which is closer to the value inferred for th e Pistol star 
(M « 4 X 10"'* Mgyr"* ; Figer et al. |l998| ), an extreme 



early- type star with similarities to rj Car. 

Since current observations and theoretical modelling 
of the optical spectrum are unable to determine the ef- 
fective temperature and the stellar radius of the primary 
without first determining Mi (cf. Hillier et al. 2001), our 



independent estimate may prove to be extremely useful 
in this regard. It will be interesting to see if the results 
from this paper are consistent with future X-ray obser- 
vations, and whether estimates of Mi from observations 
at X-ray and other wavelengths can be reconciled. Future 
X-ray grating observations should also help us to fix the 
value of the wind momentum ratio more accurately. 

As the secondary wind dominates the X-ray spec- 
trum, and its terminal velocity appears to be high 
3000 kms~^ ), we should expect to see signs of Doppler 
broadening and shifts in the line profiles. While there is 



little evidence for this in the current spectrum, other or- 
bital phases may be more favourable in this regard. The 
addition of Doppler effects has already been incorporated 
in modelled X-ray spectra (Pittard et al. in preparation), 
and should provide further information on wind velocities 
and the structure of the wind-wind collision region. 

The over-prediction of the X-ray lines in our models 
perhaps indicates that the companion has sub-solar abun- 
dances, which favours an 0-type over a WR classifica- 
tion, although we would need to perform a more detailed 
analysis to confirm this possibility. As the primary has 
slightly enhanced abundances of C and N compared to so- 
lar, this suggests that to date there has been no mass ex- 
change between the stars. Lamers et al. (199S) suggested 
that the star which dominates the UV GHRS spectrum 
is not the star which ejected the nebula since the abun- 
dances in the GHRS spectrum are not as evolved as the 
abundances in the nebula (which are indicative of CNO- 
cycle products). As the UV br ight source is probably the 
companion (Hillier et al. 2001) this indicates that it was 
the primary which ejected the nebula. There are also some 
caveats about the analysis in Lamers et al. ( |l998| ) since 
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Fig. 7. Zoomed in view of Fig. |6| showing the spectrum in the 1.8 — 2.8 keV region. The hnes are identified as either 
S or Si, and from left to right are: Si XIII triplet (1.839, 1.854, 1.865 keV), Si XIV (2.006 keV), Si XIII (2.182 keV), 
Si XIV (2.375 keV), S XV triplet (2.429, 2.447 keV), S XVI (2.622 keV). The fit to the S hnes is very good, while 
the model slightly overpredicts the emission from the Si lines, perhaps indicating that the Si abundance is slightly 
sub-solar. 



strong C lines (which Lamers et al. took to indicate nor- beyond the next periastron passage will further reveal the 
mal CNO abundances in the GHRS spectrum) also appear hidden secrets of this most enigmatic system, 
in stars known to be deficient in C (see the discussion in 
Hillier et al. 



2001). It is also interesting to note that the 



UV spectrum brightened in 1999.1 vs. 1998.2, which is 
suggestive of an eclipse of the UV source near periastron 
in 1998. In conclusion, the high energy photons (UV and 
X-ray) seem to be telling us about the companion. 



We emphasize that contrary to the vast majority of col- 
liding wind systems, our X-ray analysis of 77 Car primarily 
probes the conditions of the shocked wind of the compan- 
ion. X-ray observations of rj Car are therefore unique in 
this regard since at other wavelengths (with the possible 
exception of the far UV) the wind of the primary dom- 
inates the observed phenomena. While our analysis has 
for the first time provided a direct estimate of the wind 
parameters of the companion star, relating these to the 
stellar parameters (mass, radius, luminosity) of the com- 
panion star requires more work. It is anticipated that the 
continued multi-wavelength study of 77 Car through and 
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